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We analyze the free energy of charge and color neutral 2-flavor quark matter within the BCS approximation. 
We consider both the homogeneous gapless superconducting phase and the heterogeneous mixed phase where 
normal and BCS superconducting phases coexist. We calculate the surface tension between normal and super- 
conducting phases and use it to compare the free energies of the gapless and mixed phases. Our calculation, 
which retains only the leading order gradient contribution to the free energy, indicates that the mixed phase is 
energetically favored over an interesting range of densities of relevance to 2 flavor quark matter in neutron stars. 



PACS numbers: 25.75.Nq,26.60.+c,97.60.Jd 

O 
O 



I. INTRODUCTION 



At high baryon density, we expect on general grounds that nuclear matter will undergo a phase transition to the de-confined 
quark phase. Although the nature of this transition and its location remain unclear, theory suggests that the quark phase is likely 
to be a color superconductor. This expectation is several decades old [1, 2]. Recent work based on effective models for the 
quark-quark interaction predict a superconducting gap A ~ 100 MeV for a quark chemical potential /j ~ 400 MeV 0101 and has 
generated renewed interest in the field. The pairing or BCS (Bardeen-Cooper- Schrieffer) instability is strongest in the spin-zero, 
color and flavor antisymmetric channel, naturally resulting in pairing between quarks of different flavors. Any source of flavor 
symmetry breaking will induce a stress on the paired state since it will act to move the Fermi surfaces of the different flavors 
apart. There has been much recent interest in how this will affect the ground state properties of quark matter. In this work we 
ly-j ■ investigate the role of the large isospin breaking induced by electromagnetism in 2 flavor quark matter. 

In the absence of any flavor symmetry breaking the ground state of bulk three flavor quark matter is expected to be the Color- 
Flavor-Locked state |5]. In this phase, all nine quarks participate in pairing and the SU (3) co \ OI x SU(3)i x SU(3)r x U(1)b 
symmetry of QCD is broken down to the global diagonal 51/(3) symmetry. Chiral symmetry breaking results in light pseudo- 
Goldstone modes with the quantum numbers of the pions, kaons and eta's. This phase is neutral with respect to electric and 
color charge as it contains equal numbers of up, down and strange quarks. A finite strange quark mass breaks flavor symmetry. 
For small values of the strange quark mass (m, <C A ) the stress on the CFL state is resolved via a novel mechanism wherein 
the light pseudo-Goldstone modes condense O01- Bedaque and Schafer have shown that the stress induced by the strange 
quark mass can result in the condensation of neutral kaons in the CFL state [6]. At intermediate values of the strange quark 
mass corresponding to m 2 s //j ~ A, the situation is less clear. Recent, calculations suggest the possible existence of a color-flavor 
locked phase with non-trivial gapless quark excitations fl. 

In this article we study the case where the (effective) strange quark mass is large compared to the quark chemical potential. 
Here no strange quarks are present and the excess charge of the up and down quarks is neutralized by electrons, which in 
turn induces a splitting 8/j = (/j„ — ju<j)/2 between the up and down quark chemical potentials. If 8/j is large compared to the 
pairing energy, the pairing breaks down and the normal state is favored. At intermediate values, interesting new phases are 
possible and have been previously suggested. This includes the homogeneous gapless superconducting phase ||9|] . The gapless 
superconducting phase, also called the Sarma or breached phase| 10], was found to be an unstable maxima of the free energy 
in electron superconductors. In the quark matter context, a meta-stable gapless state was first discussed in Ref. full . In the 
condensed matter context, stable gapless phases were discussed by Liu and Wilczek fl2ll . Subsequently, Shovkovy and Huang 
showed that charge neutrality could stabilize the gapless phase in quark matter|9f]- The other possibility is the heterogeneous 
mixed phase where normal and BCS superconducting phases coexist 113111411 . Although we do not study it in this work, we 
note that there exists a small interval in 8/j/A where crystalline superconductivity could be favored llal . This crystalline phase, 
which is also called the Larkin, Ovchinnikov, Fulde and Ferrel (LOFF) phase, is characterized by a spatial variation of A on 
a scale 1 /8^. The mixed phase that we consider here shares some features with the LOFF phase: spatial variation of A and 
crystalline structure. In the mixed phase electric charge neutrality is satisfied as a global constraint- a positively charged BCS 
phase coexists and neutralizes the negatively charged normal phase. As emphasized in earlier work by Glendenning, the mixed 
phase is a generic possibility associated with first order transitions involving two chemical potentials II 1611 . This charge separation 
leads to long-range order and crystalline structure at low temperature [17]. The description and competition between the mixed 
and the gapless phases is the subject of this article. It has been shown that mixed phase is energetically favored only for small 
values of the surface tension 1 18]. In this work, we calculate the surface tension for the first time and show that it is indeed small. 
Consequently we find that the heterogeneous mixed phase is favored. 
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We will begin by deriving the free energy of color and charge neutral matter containing light (up and down) quarks and 
electrons in Sec. |I1] In concordance with earlier findings we show the existence of two stable and neutral bulk ground states: the 
homogeneous gapless state and the heterogeneous mixed phase. In Sec. IIVI we calculate the surface tension between the normal 
and superconducting phases in the mixed phase. We will then use these results to compare the free energies of the homogeneous 
and heterogeneous superconducting phases in Sec|V] Finally in Sec. IVII we conclude with a discussion of our results and explore 
some of its consequences for astrophysics of neutron stars that may contain quark matter. 



II. FREE ENERGY 



We consider a system of massless electrons, and two degenerate light quarks (up and down), with three color degrees of 
freedom at finite quark chemical potential ft = fj + hq Q + fj& Tg+fj-} Tj. We will assume that the (qq) condensation for light 
quarks is small in the quark phase and will restrict our investigation to quark masses that are small compared to the chemical 
potential. We have written the quark chemical potential in terms of the baryon chemical potential (/j = /Jg/3), the electric charge 
chemical potential /jq and the color chemical potentials /J3 and /js. The electric charge matrix and the color matrices are defined 
below. 

/ 1 n \ 1 1 / 1 \ - / 1 \ 

2=1 n 1 =Z+o T 3, T 3 =2ks= -1 , r 8 = 2V3>*= 1 . (1) 

V -3 / 6 2 \o 0/ \o o -2/ 

It is convenient to write the quark chemical potential (matrix) as 

^0 + 8^X3 +^3 r 3 +^ 8 T&, (2) 

where /jo = and 8/j = /Jg/2. The Pauli matrices T,'s act in flavor space and the Gell-Mann matrices X,'s act in color 

space. 

We model the pairing interaction by a four-quark operator with the quantum numbers of the single gluon exchange interaction. 
As Eq. l|3} shows, one gluon exchange is attractive in the color anti-symmetric channel and repulsive in the symmetric channel. 



2 (8ap8 Y 5 - 8 a5 8p 7 ) - - (8 a p8 T 5 + 8 a 8Sp Y ) 



(3) 



where the first term corresponds to the attractive color antisymmetric 3 channel. For pairing, we only consider the attractive 
channel. Further, we consider the pairing to be in the spin zero, 5-wave and flavor antisymmetric channel. So, the four-quark 
operator is taken to be 

g(¥aY5'Wp)(?^Y5'W8)(Sa7 8 p8 - 8 a5 8p 7 ), (4) 

\f = CTj7 r = iy 2 \\r*, vjF' = y T C = i^lilo- 

Using the 2SC ansatz V(?aY5 T 2¥p = /Ae a p 3 /(4g) (where blue quarks remain unpaired) for the BCS gap A in mean field approxi- 
mation, we get for the quark Lagrangian 

l +^=xj/(^ + Y0A'-w)¥ + ^(¥ay5^2Vp)(¥ T Y5W8)(Sa78pg-8 o ,58p Y ) (5) 

=tf b) («? + yew,) ¥ W - f- ~ ~ V^M^f , 
4g Z 



Hb =j"0 - 2^8 + 8^ T3 , 

¥ c ), 

- m+jo(po + 8^3 +fi%) -J5X2G2A 



*V (rg) =( \j7 y c ), 



The Dirac spinor fields \|/*), v p( rg ) represent the blue, and red-green quark fields respectively. The Pauli matrices x,-'s act in flavor 
space, and o,'s act in the red-green color space. Integrating over the fermionic variables in the partition function, the free energy 
for a constant gap field A is given by 

a 2SC (A)=^+2^l^(Vp^-^±^)Q(^T^-^p^) + y ^TrlogG(p)-^L+C, (6) 



G(p) 



■m+Yo (^0 + 8^X3 +IJ3O3 -Y5X2O2A 

Y5X2CT2A f-m- YoC"o + 8^x3 +a<303 +/j 8 ) 
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where we added the contribution to the free energy from the electrons with chemical potential /j e = —/jq- The constant C is fixed 
such that £l2Sc(A = 0) gives the free energy of the non-interacting system of electrons and quarks. 
In calculation the free energy we use the identities Jl3ll 



Trlog 



A B 
C D 



= Trlog (-BC + BDB~ l A) 



(7) 



d*p 
(2nf 



log [(p (l + i8) + 8/j + n-if - (A 2 + e 



(2^ 



{8/j + H3 + VA 2 + e± 2 ) 9 (8^ + (J3 + yjA 2 + z ± 2 ) 



- (6> + m - y/A 2 + e± 2 )9(8 J u+A'3 - + e± 2 



where K is a regularization dependent constant that does not depend on the BCS gap A. For example K = in dimensional 
regularization. In other regularization schemes such as momentum cutoff, this constant can be absorbed into the over all constant 
C defined in Eqs. (|6) and (|8j. A straightforward calculation gives 



n 2SC (A)=- 



4 



12k 2 Ag 



-2 



(Pp 
(2^ 



£ (5/j +/J3 + \/A 2 + e± 2 )9(8A( +a<3 + \/A 2 + e± 2 ) 



(8) 



v / A 2 + e± 2 )9(8 A( +A/3- VA 2 + e± 2 ) 
-p. 3 + ^A 2 + e± 2 ) 



-{8/j-/j 3 + ^/A 2 + e ± 2 )Q{5lj- 



■(Bfi-w - v / A 2 + e± 2 )9(8 A( - iU 3 - v / A 2 + e± 2 ) 



C, 



8±(p) =vT + n?±(w+w). 



The momentum integrals are evaluated with a cutoff A. The coupling g(A) is determined by requiring a BCS gap Ao = 100 MeV 
at /j = 350 MeV and /jq = = /J3 = /J8- The light quark masses were assumed to be small compared to /j and we set them equal 
to 5 MeV. 

Bulk matter is neutral with respect to electric and color charge. In what follows we investigate both the homogeneous and the 
heterogeneous phases. For the former, electric and color charge neutrality is a local condition. In this case, for a given baryon 
chemical potential /j, we require that 



2SC 



0, 



do. 



2SC 



3^3 



0: 



2SC 



(9) 



These sets of equations determine /jq, /J3, /js at any value of the gap A. Simultaneously solving the above set of equations 
together with the gap equation 3i22Sc/3A = determines the gapless states. Denoting the solutions as fiQ, £3, fig and A, the free 
energy of the gapless phase is 



^gapless (<") = ^2Sc( i ",P'e,A'3,P ! 8,A) . 



(10) 



To construct the mixed phase comprising of normal and superconducting states, electric charge neutrality is imposed as a 
global constraint. In principle we could have imposed color neutrality as a global constraint in the heterogeneous phase. How- 
ever, we expect the color Debye screening length, A-Debye-coior — 1 //A to be short and comparable to the inter-particle distance 
in strong coupling. Under these conditions, color neutrality is essentially a local constraint. For normal and superconducting 
phases to co-exist they must satisfy the Gibbs criteria: for a given /j, it requires equality of pressures at equal electric chemical 
potential. This results in the following condition 



^2SC(A = 0,/Jq) = ii2sc(A B cs,i"e) , 



(11) 



which in turn uniquely determines Jig corresponding to the mixed phase. Global charge neutrality follows as long as the normal 
and superconducting phases have opposite electric charge. The volume fraction of the superconducting phase, which we call % 
is determined by the global constraint 



where Qsf = - 



XGSF+(1-X) fiNormal =0 

3i2 2 sc(ABCS, j "e) 



(12) 



dfiQ 



2Noraial 



8n 2S c(A = 0, j u e ) 



HQ 



are the electric charge densities of the superconducting and normal phases respectively. 
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A (MeV) A ( MeV ) 

FIG. 1: In the left panel, free energy curves for the gapless (solid) and the mixed phases (dashed) as a function of A for /j = 350 MeV are 
shown. The gapless free energy is characterized by the imposition of local electric and color charge neutrality while in the mixed phase electric 
charge neutrality is imposed globally. These requirements induce chemical potentials /jq and [jg and are shown in the right panels. The filled 
dots and rectangles indicate the locations of the mixed and the gapless ground states. 

In Fig. □ we show the free energy Q.2SC as a function of A for /j = 350 MeV. The solid curve is for the case where electric 
and color charge neutrality conditions are imposed locally. The minimum seen at A ~ 40 MeV corresponds to the gapless 
superconducting state described earlier. The dashed curve shows the free energy corresponding to the mixed phase. Here we 
have imposed only local color neutrality and determined the electric chemical potential using Eq. il 11 . Accordingly, it shows the 
existence of two degenerate minima corresponding to the normal and BCS superconducting states. The variation of the electric 
and color chemical potentials are also shown in Fig. [2 We find that fii = and |/tg | <C \J1q \ < \hb\- The electric charge chemical 
potential in the mixed phase is large. For our choice of model parameters and for /j in the interval 300 — 500 MeV we find that 
the ratio Ao/8/j ~ 1.3. This is consistent with, but smaller than, the weak-coupling leading order prediction of Ao/5/j = \/2. On 
the other hand, fi% is small and proportional to Aq//j since it is only through pairing that color neutrality is upset. 

m. COULOMB AND SURFACE ENERGY IN THE MIXED PHASE 

The mixed phase free energy in Fig.^ignores finite size contributions arising due to Coulomb and surface energies associated 
with phase separation. If these effects are negligible, the results clearly indicate that the mixed phase is favored over the gapless 
phase. However, given the small difference in free energy between the gapless and mixed phases including these corrections 
becomes necessary. In Sec. II VI we calculate the surface energy in the mixed phase. For now, we treat the surface tension as a 
parameter and study how it influences the free energy of the mixed phase. 

The mixed phase can be subdivided into electrically neutral unit cells called Wigner-Seitz cells, as in the analysis of the 
inner crust of a neutron star where droplets of charged nuclear matter coexist with a negatively charged fluid of neutrons and 
electrons 1 19]. In the present context, each Wigner-Seitz cell will contain some positively charged superconducting quark matter 
and some negatively charged normal quark matter. Although at low temperature these unit cells form a Coulomb lattice, the 
interaction between adjacent cells can be neglected compared to the surface and Coulomb energy of each cell. In this Wigner- 
Seitz approximation, the surface and Coulomb energy per unit volume are fairly straightforward to calculate if the charge density 
is spatially uniform in each phase and if the surface thickness is small compared to the spatial extent of the Wigner-Seitz cell. We 
will examine both these requirements in greater detail in Sec. lIVI For now we shall assume that these requirements are satisfied. 
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In this case, the surface and Coulomb energies depend in general on the geometry and are given by 1 2' 



E> = 



d x a 
ro 



E L = 2nu em f d (x)x{AQ) 2 rl 



(13) 
(14) 



where d is the dimensionality of the structure (d = 1.2, and 3 correspond to Wigner-Seitz cells describing slab, rod and droplet 
configurations, respectively), a is the surface tension, AQ = Qsf — ^Normal is the charge density contrast between the two phases 
and oc em = 1/137 is the fine structure constant. The other factors appearing in Eqs. dl 31 and 14l > are: x, the fraction of the rarer 
phase which is equal to % for % ^ 0.5 and 1 — % for 0.5 < % ^ 1; ro, the radius of the rarer phase (radius of drops or rods and 
half-thickness of slabs); and fd(x), the geometrical factor that arises in the calculation of the Coulomb energy which can be 
written as 



fd(x) 



1 



d + 2 



'2-dx l - 2 ' d 
d-2 



(15) 



The first step in the calculation is to evaluate ro by minimizing the sum of E c and E s . The result is 



ro - 



d a 



lV3 



4%a em f d (x) (AQ) 2 _ 

We then use this value of ro in Eqs. dl 31 and dl4i to evaluate the surface and Coulomb energy cost per unit volume 



E A +E c = -Un a em d z f d {x) x z 



(16) 



(17) 



This equation allows us to include the free energy cost associated with the Coulomb field and the surface. After calculating the 
surface tension in Sec. UVI we will return to an analysis of the free energy competition between gapless and mixed phases in 
Sec.El 



IV. SURFACE TENSION 



We will now calculate the surface energy associated with the normal-superconductor interface. This will allow us to determine 
if the mixed phase is favored over the gapless state. To accomplish this we will need to allow for spatial variations of A and 
calculate the gradient contribution to the free energy. Near the critical temperature one can employ Landau-Ginsburg theory 
to obtain the leading gradient contribution to the free energy. In this case, the gradient expansion (in powers of |V A(x)| 2 ) 
is well defined. Bailin and Love have derived an expression for the gradient contribution to the free energy for a relativistic 
superconductor near the critical temperature 01 . In the context of den se q uark matter, a discussion of the Landau-Ginsburg 
theory for two and three flavor quark matter can be found in Refs. I22T l23ll . For temperatures that are small compared to the 
critical temperature, the gradient expansion (in powers of | V A(x) | 2 ) has poor convergence properties and one must in principle 
include gradients to all orders. However, if the spatial variations of A occur on a length scale that is numerically large compared to 
A -1 , it is a good approximation to retain only the leading term. We shall assume that the thickness of the normal-superconductor 
interface is larger than A -1 and retain only the leading contribution from the gradient term. In this case, the A dependent 
contribution to the free energy of the 2SC phase, including the gradient term, is given by 

J (A) 2SC =02SC + K§c|V A(*)| 2 + 0(|V A(x)| 4 /A 4 ) . (18) 

The coefficient of the gradient term (K2sc) for the 2SC phase, which is related to the Meissner mass of the gluons, at zero 
temperature and for the case where 8/j = has been calculated in earlier work by Rischke 12411 . Recently, Huang and Shovkovy 
have investigated how 8/j affects the Meissner masses (which as discussed earlier, in turn determines K2sc) ES- They find that 
the Meissner mass of the eight gluon (assuming that the condensate aligns in the 3 direction) remains unaffected by 8/j in the 
BCS phase i.e, for 8/j < Ao. On the other hand, they find the Meissner masses of gluons with color 4,5,6,7 (which are degenerate) 
decrease with increasing &/j and become negative when 8/j = Ao/ v2- Surprisingly, these gluon masses become negative even 
in the BCS phase, i.e, for Aq/^/2 < 5/u < Ao . It is presently not clear how this instability is resolved. In this work, where we 
have employed an Nambu-Jona-Lasino (NJL) model description of the interaction between quarks, the gluons are introduced as 
fictitious gauge fields solely for the purpose of determining the coefficient K2sc in the effective theory for the real field A I25ll . 
By minimal coupling between the field A and the gluons, the relation between K2sc and the Meissner mass (m 2 ~ CC 2 /j 2 ) of the 

gluons is obtained by matching to be K 2 ^ c m^/oc^A 2 ,, where a s is the strong coupling constant. 
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FIG. 2: The profile of A(x) at the interface for /j = 400 MeV. The solid curve shows the variation of the dimensionless quantity §(x) = A(x)/Ao 
as function of the x = xAq. The dotted-line corresponds to the region where the leading order gradient energy is less than zero when 8/j 
corrections from Ref. |21] are included. The dashed curve shows the variation of VA(x)/A(x) 2 . 



(2) 

In our calculation, the coefficient is determined by matching the mass of the fourth and eight gluon in the effective theory 

(21 

to the microscopic calculation. We note that Kjer could depend on A(x), since the low energy effective theory for the real field 
A(x) does not have a well defined expansion parameter. Hence, in general 

( 2) y (A(X)^ ^ 

n =o,~ V A 

where a n are constants that remain to be determined by matching conditions. We make the simplifying assumption made in 
Ref. 1 24] and retain only the leading (ao) and next to leading {a\) order terms. In this case, using the Meissner masses calculated 
in Ref. 12 ill and for Aq > 8/j we obtain 



A 2 / 2 V Al> A 2 



K2SC "37r2A2 



(20) 



where p = /j+/jq/6+/j^ is the common chemical potential of the quarks participating in pairing and S/j =/jq/2. From Eq.l20lwe 
see that for Aq/S/j < \fl, the coefficient ao < while a i remains positive in the mixed phase where Ao /8/j > 1 . In weak coupling 
and to leading order in the Ao and 8/j, the mixed phase occurs precisely at Ao/8/j = \fl. In our calculation of the mixed phase free 
energy in Sec.|II|we found that Ao/8/j ~ 1.3 for /j in the interval 300 — 500 MeV. This result, which implies that Ao/8/j < is 
not generic to the mixed phase. For other values of the coupling we indeed find that Ao/8/j > \/2 and both ao and a i are positive. 
Nonetheless, even when ao < 0, for A(x)/Ao ~ 1, the a\ term ensures that the gradient contribution is positive and is therefore 
crucial for the stability of the BCS phase component of the mixed phase. In the vicinity of the BCS phase, i.e, when A(x) /Ao ~ 1 , 

the gradient contribution is independent of 8/j and = /j 2 /67I 2 Aq. We note that had we retained only the ao contribution to 

(2) 

K 2SC and matched to the mass of the eight-gluon (arguing that in our NJL model treatment only the generators T3 and 7g are 
relevant) we would still have obtained K^ c — /j 2 /67I 2 Aq. A non-zero 8/j in general acts to reduce the energy cost associated with 



the gradient term when A(x) /Aq decreases from unity. Surprisingly, when A(x) /Aq = J 2(2(8 J u 2 /A 2 1 ) — 1 )/ (4(8/j 2 / Aq) — 1 ), the 
gradient energy goes to zero. This does not necessarily mean that the interface is unstable since we have not investigated the 
higher order gradient terms. However, this trend suggests that we can obtain an upper bound on the surface tension by using the 
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(largest) value of K^gr = /U 2 /6ji 2 Aq across the whole interface. We note in passing that for temperatures near 7c, where Tc is the 
BCS critical temperature, Bailin and Love [2] find that 



,(2) _ 7C(3) 



V 2SC 



24ti 4 a 2 A 2 , 



(21) 



(2) 

where a ~ 0.57 = kgTc /Ao, C(3) — 1-202. Interestingly, the results for Kio C in the vicinity of the BCS state at T = and T ^Tq 
are both parametrically and numerically similar. 

The contribution to the free energy which is independent of spatial derivatives is given by i22sc(A) and was derived earlier in 
Eq. In general, both the gradient term and i22sc WQ functions of all three chemical potentials, namely: /j, /jg and /jg. While 

we retain this full dependence in £l2SC we use the approximate form for K 2SC valid when A(x) ~ Ao. To construct the spatial 
profile of A between the normal and superconducting phases we use the equation of motion. In one spatial dimension this is 
given by 



d 2 A(x) dA(x) 2 3k 2S c 
K2SC — — — I- ( — ; — ) 



2SC 



dx 2 



dx 



dA 



dA 



(22) 



We solve this ordinary differential equation subject to the boundary condition that A(x = 0) = and A(x) = Ao(u) at large x, 
where Ao(p) is the gap in the BCS phase. The results are shown in Fig.|5J where we calculate the interface for various values 
of the baryon chemical potential /j = 400 MeV. We find that the profiles for the interface for other values of /j in the range 
350 — 450 MeV show nearly identical behavior when plotted in terms of the reduced dimensionless units §(x) — A(x)/Ao(/j) 
and x — xAq(jj), The typical length scale for the variation of A is parametrically similar to but numerically larger than Aq . We 
find that a rough measure of the thickness f ~ (1.5 — 2)Aq . As mentioned earlier, due to the lack of a well defined gradient 
expansion for effective theory at T = 0, the validity of the leading order (in the gradient expansion) results requires that the 
spatial variation of A(x) be weak and satisfy the condition VA/Aq -C 1. From Fig. |2] we see that this condition is satisfied to a 
fair degree (see dashed lines). Using these profiles, we have computed the surface tension by integrating the free energy function 
given in Eq. Jl 81 . The results are presented in the table below. 



fi (MeV) 


O (MeV/fm 2 ) 


H (MeV) 


O (MeV/fm 2 ) 


fi (MeV) 


O (MeV/fm 2 ) 


350 


1.9 


375 


2.6 


400 


3.6 


425 


4.9 


450 


6.8 


475 


9.4 



(2*) 

Our estimation of the coefficient of the kinetic term, which we call K 2 § c ignores several possible corrections. Including those 
due to the truncation of the series in Eq. \19\ and the neglect of 8/j. As discussed earlier the effect of 8/j is to decrease the 
gradient contribution, consequently we can expect these corrections to lower our estimate of the surface tension. In Fig. |2] the 
vertical dotted-line indicates the region where the gradient energy is less than zero when 8/j corrections are included. In this 
case, the profile would resemble the dotted-line rather than the solid line and corresponding surface energy is greatly reduced. 

(2) 

We have also examined how an increase in kv c would affect our calculation of the surface tension. We find that these changes 
are modest. The surface tension a increased by 40% when K2sc was increased by a factor of two. 



V. PHASE DIAGRAM 



Using the results of the previous sections, we can study the free energy competition between the gapless and the mixed phases. 
The free energy of the mixed phase, including Coulomb and surface contributions is 

Cl mixe d(lj)= discing) +E S + E C , (23) 

where the surface and the Coulomb contributions were defined in Eq. dl7> and Jig is defined through Eq. il It . The left panel 
of the Fig. [5] shows the free energy difference Af2 = fi ga pi eS s — n m i xe d for various values of the surface tension in the mixed 
phase. The filled circles correspond to the values of the surface tension computed in Sec. II VI We have used this information in 
constructing the "phase diagram" as a function of /j and o, and is shown in the right panel of the figure. As expected, the gapless 
homogeneous phase is favored when a is large and the mixed phase is favored for small a. The surface tension computed in the 
previous section is also shown (filled circles). Their low values at small /j indicate that the mixed phase is favored at low density. 
At first, with increasing density the free energy difference (A£l) increases, resulting in a robust mixed phase. For /j> 450 MeV 
the free energy difference decreases, indicating a possible first-order transition to the gapless phase at high density. The dashed 

(2) 

dashed-curve on the right panel shows the surface tension for the case where we artificially increased kW c by a factor of two. In 




FIG. 3: Left panel: Influence of the surface energy on the total free energy difference between the gapless and the mixed phases Afl = 
f^gapless ~ ^mixed f° r various values of the surface tension. The filled circles correspond to values of o(fi) shown in the right panel. Right 
panel: The phase diagram showing the favored phase as a function of the surface tension and baryon chemical potential. The filled circles 
correspond to the results obtained from calculation of the surface tension outlined in Sec. lIVI The dashed line shows the surface tension when 

(2) 

K2 SC was increased by a factor of two. 



this case the trends are similar, the mixed phase continues to be favored at low density and a first order transition to the gapless 

(2) 

phase occurs at /j ~ 500 MeV. To eliminate the mixed phase l&ar has to be at least on order of magnitude larger than the one 
used in this work. 

We now examine if the assumptions made earlier in constructing the mixed phase are satisfied for the values of the surface 
tension calculated in Sec. IIVI The Debye screening length in the normal and superconducting phases are easily computed by 
using the relation A,^ 2 = — 4% a em 3 2 ^2Sc/3j"q- For the individual phases to have uniform charge distributions, their spatial 
extents in the mixed phase must be small compared to the Debye screening length. The table below provides a comparison 
between these length scales at various values of the baryon chemical potential. The results are shown for the mixed phase with 
d = 3, corresponding to a spherical geometry for the Wigner-Seitz. In the table, Abcs, X, Rd and R ws correspond to the gap in 
the BCS phase, the volume fraction of the BCS phase, the droplet radius and the radius of the Wigner-Seitz cell respectively. 
Further, we have also implicitly assumed that the interface thickness is small compared to Rd and R ws . The thickness parameter 
f, defined to be the distance over which the gap roughly decreases by a factor of e ~ 2.718, is also given in the table. 



H (MeV) 


Abcs (MeV) 


X 


Rd (fin) 


Rws (fm) 


A,g (fm) 


^normal (fm) 


t (fm) 


350 


74 


0.4 


6 


8 


7 


5 


4 


400 


106 


0.7 


4 


6 


6 


4.5 


3 


450 


140 


0.8 


3 


5 


5 


4 


2 


500 


174 


0.9 


3 


5 


5 


4 


2 



The length scales shown in the above table are just barely compatible with our earlier assumptions. In particular, the Debye 
screening length and the spatial extents are of similar size, the charge distribution in the mixed phase could differ from the 
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simple profiles considered here. We also note that the thickness is not small compared to the typical size of the droplet or the 
Wigner-Seitz cell. We will not consider these finite size corrections in this work. These corrections generically tend to decrease 
the Coulomb energy, increase the surface energy and penalize the mixed phase when the volume fraction % is near zero or unity 
ll26ll27ll . These finite size effects in the mixed phase clearly warrant further study. 




FIG. 4: A schematic phase diagram at mixed fj. The BCS gap Ao (or equivalently the effective four-fermion coupling constant) is plotted 
along the x-axis and d/j = jjq/2 is the y-axis. The line 8/j = Arj/\/2 separates the regions where the BCS state (lower-right) and the normal 
states (upper-left) are favored. The charge-neutral ground state follows the trajectory shown by the thick lines with arrows. The wedged region 
between the lines corresponding to 8/j = Ao/\/2 and 8[j = 0.75Ao indicates the region were the LOFF phase may be favored| 15]. 

The results of this work combined with the finding of earlier investigations in Refs. fl^.l28ll29ll32il of the charge neutral two 
flavor quark matter suggest a schematic phase diagram of the form shown in Fig.|4] The figure depicts the phase structure at 
fixed /j . The BCS gap Ao or equivalently the effective four-fermion coupling increases along the x-axis and 8/j — hq/2 increases 
along the y-axis. The line S/j = A/V2 separates the regions where the BCS state (lower-right) and the normal states (upper-left) 
are favored (As discussed earlier, this is true only in weak coupling. In strong coupling, 8/j lies in the vicinity of Ao/y/2). The 
charge-neutral ground state follows the trajectory shown by the thick lines with arrows. In each phase, regions above this line 
correspond to net positive charge and regions below to net negative charge. For small coupling (Ao), the homogeneous normal 
phase is preferred, and is represented by the line from point A to B. At intermediate coupling, the mixed phase containing an 
admixture of negatively charged normal phase coexisting with positively charged BCS phase is favored. This is shown by the 
line from point C to D. At large coupling, the homogenous BCS or 2SC phase prevails and is depicted by the line from point D 
to E. The wedged region between the lines corresponding to 8/j = A/V2 and 8/j = 0.75A indicates the region were the LOFF 
phase may be favored 1 15]. It remains to be seen how this phase could exist as a charge-neutral state. If the charge density in the 
LOFF phase resembles the normal phase, the figure supports the existence of a region between point B and C where the LOFF 
phase could be stable and electrically neutral. It also points to the possibility of a mixed phase of LOFF and BCS. Investigation 
of the LOFF phase in the charge neutral phase diagram is clearly warranted but is beyond the scope of this article. 

VI. DISCUSSION 

Based on our calculation of the surface tension between normal and superconducting 2-flavor quark matter we have shown 
that a heterogeneous mixed phase is likely to be the ground state of quark matter at low density. The stress induced on the bulk 
system by the requirement of charge neutrality is resolved by phase separation. The other possibility, which is the homogeneous 
gapless phase is energetically disfavored at low density. As depicted in Fig. |4] at small coupling we can expect the charge 
neutral ground state to be be in the normal phase. At intermediate coupling, the heterogeneous mixed phase is favored and large 
coupling a uniform BCS like 2SC phase is the ground state. 

While we expect our qualitative results to be fairly robust, there are three important caveats to our quantitative findings that 
we shall now discuss. First of these caveats is related to the earlier observation that the surface thickness t ~ 1 — 3 fm, the typical 
physical sizes of the charged regions ~= 2 — 6 fm and the Debye screening length Xd — 5 — 8 fm are not well separated. 
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This is likely to alter the simple structure of the mixed phase considered in this work. The charge, baryon density and the gap 
will vary smoothly over most of the Wigner-Seitz cell. These finite size corrections will need to be examined before one can 
draw firm conclusions regarding the phase competition (for a discussion of these effects in the context of a nuclear-kaon mixed 
phase see Ref. iS). Second, we have assumed that the weak coupling BCS theory holds even when A//j ~ 1/3. In the BCS 
limit A -C /J, neither the gapless nor the mixed phase could exist because for 2-flavor quark matter the electric charge chemical 
potential is proportional to the baryon chemical potential. In the strong coupling regime it is difficult to know a priori how 
other contributions to the free energy would alter this competition. Finally, we have only included the leading order gradient 
contribution in the effective theory for A. Although our profiles indicate that spatial variation of A(x) is mild (From Fig. [2] 
|VA(x)| 2 /Aq < 1), it is a priori not possible to assess how the higher order spatial derivatives of A(x) in the free energy would 

(2) 

change our estimate of the surface tension. As could corrections to the coefficient K 2SC arising from strong coupling and the 
lack of well defined expansion parameter for the effective theory at T = 0. Additionally, corrections due to electric and color 
charge chemical potentials are important [21 ]. Interestingly, the correction to the gradient term from the electric charge chemical 

(2) 

potentials tends to decrease the surface energy. We have shown that drastic changes to K 2SC (an increase by a factor of 2) resulted 
in a surface tension that continued to favor the mixed phase at low density. In this case, with increasing density, we found that 
the gapless phase becomes favored. The location of this first order transition between the mixed phase and the gapless phase is 
model dependent since it is sensitive to the density dependence of the ratio /j/Aq. 

An important consequence of 8/j corrections to the kinetic term concerns stability of the gapless phase. First we note that 
we can expect the matching between the kinetic coefficients for the effective theory for A(x) and the masses of gluons in the 
microscopic theory to be valid in the gapless phase. Since 8/j > Ao in the gapless phase this implies an instability of the spatially 
uniform gapless phase (the leading order gradient energy is negative). How this instability is resolved is not clear at this time. It is 
however clear that some type of heterogeneity must occur. Likely candidate states include the LOFF phase and the mixed phase. 
Pending a detailed investigation of the charge neutral LOFF state, our results suggest that this instability may be resolved by the 
formation of a mixed phase. In our NJL model study, the mixed phase is stable with respect to small gradient perturbations for 
all Ao/8/j > 1. In weak coupling, since Aq/8/j = \/2 in the mixed phase, the chromo-magnetic instability discussed in Ref. [21] 
for gluons 4-7 in QCD does not occur in the mixed phase. In strong coupling, whether this instability persists and if so, the 
nature of its resolution is not yet clear. 

A heterogeneous quark matter ground state has important consequences for neutron star physics - were quark matter to exist 
inside their cores. The free energy difference between the gapless and mixed phase is small. Consequently, its effects on neutron 
star structure will be small. On the other hand, we can expect important changes to the transport properties induced by the 
heterogeneity of the mixed phase. In this phase, transport of the low energy neutrinos, photons, etc., will be dominated by 
coherent scattering induced by the presence of superconducting "bubbles" of quark matter with typical dimension of ro ~ 5 fm 
(for a discussion of how heterogeneity affects low energy neutrino transport see Ref. |30[). Further, the existence of regions 
with normal quark matter will affect the neutrino emissivities - since they are typically large in the normal phase. Finally, the 
crystalline structure of the mixed phase could be relevant to the modeling of glitches observed in some pulsars. 

We note that there are other mixed phases that could play a role in neutron stars. These include the mixed phase between quark 
matter and nuclear matter| 16]; between color-flavor-locked matter and nuclear matter|31], or superconducting quark matter and 
normal quark matter|32]. All of these share several common features: the coexistence of positively and negatively charged 
phases separated by an interface; crystalline structure at low temperature; and coincidentally typical droplet sizes of the order 5 
fm. However, in these aforementioned examples it is still unclear if the mixed phase is favored when surface and Coulomb costs 
are included. This is because the surface tension is poorly known. Arguments based on naive dimensional analysis indicate that 
the surface energy cost may indeed be too large in these systems 13 lL 13 3tl . The exception could be the mixed phases of normal 
and superconducting quark matter considered in Ref. 13211 - In this case, the surface tension is computable and the physical 
considerations are very similar to those studied in this work. 

Although quantitative aspects of the analysis presented in this work pertains to asymmetrical relativistic superconductors, 
our finding provide some qualitative insight into the behavior of non-relativistic superfluid that can be realized in ultra cold 
fermionic-atom traps being developed in the laboratory. Although there is no Coulomb cost associated with the mixed phase 
in these systems, a surface energy cost must be overcome when the surface to volume ratio is not negligible. The methods 
described in Sec. lIVI applv in general to asymmetrical fermion systems and its relevance to laboratory experiments is currently 
under investigation. 
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